IPHT-T09/196 
CERN-PH-TH-2009-230 



Geometrical interpretation of the topological recursion, and 
integrable string theories. 



Bertrand EynarJ^, Nicolas Orantii^ 

Institut de Physique Theorique, 
CEA, IPhT, F-91191 Gif-sur-Yvette, France, 
CNRS, URA 2306, F-91191 Gif-sur-Yvette, France. 

Theory division, CERN 
CH-1211 Geneva 23, Switzerland. 



Abstract: Symplectic invariants introduced in [23] can be computed for an arbitrary 
spectral curve. For some examples of spectral curves, those invariants can solve loop 
equations of matrix integrals, and many problems of enumerative geometry like maps, 
partitions, Hurwitz numbers, intersection numbers, Gromov-Witten invariants... The 
problem is thus to understand what they count, or in other words, given a spectral 
curve, construct an enumerative geometry problem. This is what we do in a 
semi-heuristic approach in this article. Starting from a spectral curve, i.e. an 
integrable system, we use its fiat connection and fiat coordinates, to define a family of 
worldsheets, whose enumeration is indeed solved by the topological recursion and 

symplectic invariants. In other words, for any spectral curve, we construct a 
corresponding string theory, whose target space is a submanifold of the Jacobian. 
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1 Introduction 

Topological String Theories aim at addressing the question of "counting" how many 
Riemann surfaces (worldsheets) with given boundary conditions, can be embedded 
into a given target space. Witten suggested [H] an underlying string field theory, in 
which worldsheets are obtained by gluing some basic building blocks. For example in 
Teichmiiller theory, building blocks are "pairs of pants". In Kontsevich's approach, 
building blocks are cylinders glued along a ribbon graph, and this idea has then given 
many variants. 

Recently, it was suggested by BKMP [TU], that Gromov- Witten amplitudes of the 
type A topological strings in a toric CY 3-fold target space X, coincide with the "sym- 
plectic invariants" (introduced in [23]) of the spectral curve of the mirror X of the 
target space X: 



The main interest of that conjecture, is that the right hand side, i.e. the symplectic 
invariants, is much easier to compute than the left hand side for given genus. 

Here, we shall study this claim, and try to understand its geometric meaning. 

In fact, we shall work backwards, and starting from a spectral curve S, we shall 
try to construct a "string theory" whose partition function is given by the symplectic 
invariants Fg's. 

The basic idea, is that out of a spectral curve, we can construct an integrable 
system, and in particular fiat connections, and a system of action-angle variables. For 
every initial condition, the angle variables of an integrable system, follow a uniform 
linear motion in time, which means that they generate a 1-dimensional manifold in the 
phase space, and, because of conformal invariance, time is a complex variable, it means 
a 1-dimensional complex manifold, i.e. a Riemann surface embedded in the phase space 
(the phase space of action angle variables has a toric symmetry). Moreover, since the 
motion is uniform, the angle coordinate is a flat coordinate almost everywhere on each 
such Riemann surface, and this gives a natural foliation on all surfaces. These surfaces 



BKMP conjecture: 



GW,(X) = F,{S^). 



(1-1) 
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can thus be cut into "propagators" and cylinders in a unique way. This immediately 
implies that the generating functions which count such Riemann surfaces of a given 
topology, do satisfy the topological recursions of [23] , and thus they are the symplectic 
invariants. 

Outline 

In section 2, we recall the definition of the symplectic invariants as well as some of 
their properties. We also recall that they have a diagrammatic representation, which 
resembles very strongly what could be expected for a string field theory. 

In section 3, we consider an arbitrary integrable system, we recall the notion of 
"action-angle" flat coordinates, and how action angle coordinates can be used to gen- 
erate worldsheets of a string theory. In other words, we define an adhoc string theory 
attached to an integrable system. 

In section 4, we define moduli spaces of worldsheets of given topologies and brane 
boundary conditions, and we show that worldsheets can be decomposed into propaga- 
tors and cylinders. This induces a decomposition of moduli spaces of worldsheets into 
cells labeled by graphs. 

In section 5, we translate the decomposition of moduli spaces in terms of string 
amplitudes, in local patch coordinates. A consequence is that, after Laplace transforms, 
string amplitudes obey the topological recursion of [23]. In particular, closed string 
amplitudes of genus g are the symplectic invariants Fg. 

In section 6, we rewrite amplitudes in terms of intrinsic geometry of the spectral 
curve. That allows to identify the spectral curve with the disc amplitude, and the 
Bergman kernel with the cylinder amplitude. 

In section 7, we show how expanding the generating functions in terms of other 
formal variables, through Lagrange inversion formula, can give to many combinato- 
rial identities. This generalizes Cut and Join equations and ELSV formulae. Indeed, 
expansions near branchpoints of the spectral curve can always be written in terms of 
intersection numbers, and expansions near singularities of the spectral curve can be 
written in terms of winding numbers. 

In section 8 we discuss the case of toric CY target spaces leading to a geometrical 
interpretation of BKMP conjecture. 

Section 9 is the conclusion. 
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2 Symplectic invariants of a spectral curve 

Let iS be a spectral curve, i.e. it is the data of a compact Riemann surface C of genu^ 
Q, and two analytical functions x and y on C, or on some open domain of C: 

S={C,x,y). (2-1) 

Typically, in topological strings, C is an algebraic curve of equation iJ(x+,a;„) = 
where H is some polynomial, and x = In (a;_|_), y = \n (x_). 

Notice that (C, x) is a Hurwitz space, i.e. the data of a compact Riemann surface 
C together with a projection x : C — CP^, which realizes C as a branched covering of 
CP^, with a branching structure given by the zeroes and poles of the differential dx: 
the branchpoints are the zeroes and poles of dx, and the monodromies are the orders 
of the zeroes and poles of dx. 

From now on, we assume that dx is meromorphic, and all zeroes of dx are simple 
(this is the case for spectral curves of topological strings, and also for spectral curves 
of matrix models counting discrete surfaces). 

• Bergman kernel 

Then we define a Bergman kernel B{zi,Z2) on C, i.e. a 2nd kind meromorphic 
symmetric 2-form on C having a double pole with vanishing residue at zi = Z2 and no 
other pole. It is normalized by requiring that near Zi = Z2, in any parametrization $,{z) 
it behaves like: 

B{zuZ2) = B{z2,z,) ~ ^iifil^^ + regular. (2-2) 

' ^^2 (^(2:1) -^{Z2)y 



zi- 



As defined here, the Bergman kernel is not unique, one may add to it any combination 
of holomorphic forms, i.e. differential forms without poles. 

For a given symplectic basis of noncontractible cycles on C, i.e. 2q cycles {Ai, B^}^^-^^, 
satisfying: 

one can define a basis of holomorphic forms {dui}f^i normalized by 

(p duj{z) = 6ij and (p duj{z) = Tij (2-3) 

where r is the Riemann matrix of periods (see for instance [271 [2H])- One can then 
parameterize the holomorphic deformations of the Bergman kernel with a symmetric 



^the genus q of the spectral curve has nothing to do with the genus g of worldsheets studied further. 
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matrix k of size g x g, and one may consider the Bergman kernel shifted by a combi- 
nation of holomorphic forms, as a new admissible Bergman kernel: 

9 

B{zi, Z2) B{zi, Z2) + 2m ^ dui{zi) duj{z2). (2-4) 

A choice of k is more or less equivalent to a choice of a symplectic basis of cycles 
Ai{K), Bi{K,), z = 1, . . . , g on which the Bergman kernel is normalized by 



/ 

J zieAi(K) 



B{zi,Z2) = 0. 

From now on, let us assume that we have chosen a Bergman kernel, or in other words 
a matrix k, or in other words a basis of cycles. 

• Branchpoints and conjugated points 

The branchpoints are the points with a "vertical tangent", i.e. the zeroes of dx: 

dx{ai) = 0. (2-5) 

We assume that all branchpoints are regular, i.e. they are simple zeroes of dx, and 
they are not zeroes of dy nor poles of y. This means that near a branchpoint Zi, the 
curve behaves like a square- root: 

y{z) ~ y{ai) + Ci^ x{z) - x{ai) ^ ... (2-6) 

This also means that in a small vicinity of Oj, there is a unique point z ^ z\\y the same 
vicinity of such that: 

x{z) = x{z). (2-7) 

z is called the conjugated point of z. It is defined only locally near branchpoints, and 
it is not necessarily defined globalljfl. 

• Recursion kernel 

We define the recursion kernel: 

/a.:..)=/a..fl = ^^j^i||g^. (2-8) 

It is a 1-form in Zi, defined globally on zi G C. It is the inverse of a 1-form in z, defined 
only locally near branchpoints. 



■*A notable exception is the case of hyperelliptical surfaces, where z ^ z is the hyperelhptical 
involution and is defined globally. 
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Figure 1: Branchpoints are points with a vertical tangent. Near a branch point 
there are two branches coming together. If z is on one branch, we call z, the point on 
the other branch with the same x projection x{z) = x{z). Notice that z is not globally 
defined, it is defined only locally near branchpoints. For example, if z moves from ai 
to 02, the analytic continuation of z would take the wrong branch. 



It has the property that it has a simple pole at branch points, and near z — > it 
behaves like: 



• Symplectic invariants 

Then, following [23], we define a sequence of symmetric meromorphic n-forms, called 
Un^ for every n and g integers, by the following recursion (often called "topological 
recursion" ) : 



K{z^,z) 




+ regular. 





(2-10) 
(2-11) 






,(9-h) 
l+n-\I 



(2, J\I) 



h=0 ICJ 



(2 - 12) 



where J is a collective notation J 



{^2, • • • , Zn+i}, and Yl' means that we exclude the 



terms (/i, /) = (0, 0), ((?, J). 
We also define ii g > 2: 




(2-13) 
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where $ is any analytical function defined in the vicinity of branchpoints such that 
(i$ = ydx. There are also definitions for Fq and Fi, and we refer the reader to [23] for 
those two cases. Fq is often called the prepotential. 

All the a;i^^'swith2-2^-n < are called stable, and the others are called unstable. 
The only unstable ones are Fq, ujf\u)^\ Fi. 

Although the definition doesn't look symmetric, every ujn'^ is a symmetric n-form. 
Stable c^i^^'s have poles only at branchpoints, of order at most Qg — A + 2n, and with 
vanishing residues. 

2.1 Some properties of symplectic invariants 
• Rescaling 

Under a rescaling y Xy, we have (if 2 — 2(7 — n 7^ 0): 

cu^5) ^ A2-23-"^te), (2-14) 

and in particular (for g ^ 1): 

Fg^X'-'^F,. (2-15) 
In particular, this implies that Fg is invariant under the parity transformation y — > — y. 

• Symplectic invariance 

If two spectral curves S = (C, x, y) and S = (C, x, y) are such that there is an 
analytical bijection from C ^ C which conserves the symplectic form in C x C: 

dx Ady = dx A dy, (2-16) 

then we have (for g > 2): 

Fg{S) = Fg{S). (2-17) 
This is why we call the Fg{S) the symplectic invariant of degree 2 — 2g. 

• Other properties 

There are many other properties, for instance concerning modularity, infinitesimal 
variations of spectral curve, singular limits, and integrability. For example, for any 
spectral curve, the F^'s satisfy holomorphic anomaly equations [22] . 

Also, it turns out that the F^'s allow to construct a Tau-function, and an integrable 
system associated to S. 

All those properties can be found in [23l [26] . 
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2.2 String Field theory: diagrammatical rules 

Let us represent pictorially every Un^ as a surface of genus g with n boundaries labeled 
by zi, ...,Zn. 




Let us represent our two kernels as elementary pieces used to build such surfaces: 
• The "propagator" 




(2-18) 



• The "cylinder": the two point function (2:1,^2) 




(2-20) 

where one integrates over the intermediate variable z. This recursion is said to be 
topological since the surfaces generated in the right hand side have strictly higher Euler 
characteristic than the one of the left hand side, and thus this recursion terminates after 
a finite number of steps equal to minus the Euler characteristics. 

In other words, the topological recursion tells us that every surface enumerated by 
the Un^^s, can be decomposed into 2g — 2 + n propagators and n + g — 1 cylinders, 
more or less in a unique way. 

Example: 
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For Ui)^\ our recursion formula gives in two iterations: 





(2-21) 

which is merely the pictorial representation of the following residue formula: 
^^\zi, 22) = 2 V Res Res K{zi, z)K{z, z') \b{z', Z2)B{z, z!) + B{z, Z2)B{z', z! 



That diagrammatic representation of the topological recursion, resembles strongly 
what one could expect to be "string field theory" diagrammatic rules, i.e. how to 
compute string theory amplitudes by gluing surfaces. 

In fact it was conjectured by BKMP [10] that Gromov-Witten amplitudes, i.e. 
topological string theory amplitudes, which enumerate surfaces embedded in a certain 
target space X, should be equal to the symplectic invariants of the spectral curve which 
is the singular locus in the mirror of the target space X. Here we are going to try to 
explain heuristically why the topological recursion indeed counts some string theory 
amplitude by giving a meaning to the graphical representation of the present section. 

But we are going to proceed backwards, i.e. given a spectral curve and its symplectic 
invariants, we are going to construct a corresponding effective string theory. 

2.3 Flavor of string theory 

Just in order to give some intuition of what we do in the next section, we recall in a 
very sketchy way, the basic ideas underlying string theory. 

Consider a "space-time" called target space X, which is a complex manifold of 
fixed dimension D. 

A "closed string" is an embedding in X of a circle: |x(cr) , a G [0, 27r]| C X with 

X{0) = X(27t), and an open string is an embedding of a segment, i.e. we don't assume 
X{0) = X{27i). From now on, let us consider only closed strings. As the string moves 
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with time T in the target space, the history of the string sweeps a surface, called world 
sheet in the target space: 

X{a,T)eX (2-22) 

where a G [0, 27r] is the coordinate along the string, and T G M+ is the time coordinate, 
and X = {Xi, . . . , X^) is a point of X, in a local coordinate system. 

String theory amplitudes are obtained by "counting" how many histories can relate 
an initial state to a final state, i.e. enumerating surfaces having given boundaries. 
Surfaces should also be counted with a weight, typically the exponential of some action 
(for instance Nambu-Gotto's action depends only on the total curvature and area of the 
surface spanned in the target space). In other words, we have to perform a functional 
integral over all coordinates Xj((T, T)§|: 

r ^ 

(boundaries) ":=" / JJ P[Xi((T, T)] e-^="°"[^('^'^)l. (2-23) 

1=1 

The "functional measures" V[Xi{a,T)] can be quite complicated depending on the 
geometry of X. When X is a submanifold of C"^^^', the functional integrals can be im- 
plemented in terms of standard functional measures in C^+^', by Lagrange multipliers 
which enforce D' relationships between the components Xi, i = 1, . . . , D + D'. This 
procedure leads to a so-called cx-model description of this string theory. This is just 
one representation of that theory among others but it has the advantage of being suffi- 
ciently well known, for some classes of target spaces, to be mapped to some integrable 
system (see [T7] for example). The boundary conditions can also be implemented by 
Lagrange multipliers. 

Now we shall assume that the theory has a conformal invariance property, i.e. that 
the functional measures and the action, are invariant under conformal reparametriza- 
tions of the worldsheet. In other words, we want to count only once worldsheets which 
are conformally equivalent. Choosing one representant per equivalence class is often 
realized by "gauge fixing", i.e. by introducing "ghosts", but we shall not need it in 
this article. 

We shall only say that X{a, T) = X(cr + iT) is an analytical function of a complex 
variable t = a + iT, and conformal reparametrizations, are obtained by changing t to 
an analytical function of t. In other words, the worldsheet is a 1 dimensional-complex 
motion X(t) in X with a complex time t. Gauge fixing means choosing a complex 

^Following our description, these amplitudes are scattering of closed strings. However, they are 
referred to as open amplitudes in the topological string literature since they enumerate surfaces with 
boundaries, i.e. open surfaces. This name can also be seen as originating from the open-closed duality 
obtained by considering a as the time instead of T. The worldsheet is then spanned by open strings 
ending on some manifolds called Branes. 
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coordinate t on each worldsheet. One of the main points in our article, is to find a 
"canonical" choice of time coordinate. 

Then, we shall make another assumption, which is that our action is integrable. 
There are many definitions of integrability. At the classical level, one of them can 
be phrased like: any motion which extremizes the action (and then called "classical 
motion"), has as many conserved quantities as half the dimension of the phase space, 
and which thus implement the same number of commuting hamiltonian fiows, and 
thus imply a toric symmetry. Another one, more convenient for us, is that, for every 
classical motion, there exists a suitable change of variables, mapping the coordinates 
Xi{t), Xi{t) to the so-called "action and angle" variables such that after the change of 
variables, the motion is linear at constant speed in a mult i- dimensional torus (see [7| 
for an introduction to integrable systems). This also provides a natural torus action. 

A consequence of integrability and a torus action, is that there is a localization 
formula, and the functional integral above can be reduced to a sum over only world- 
sheets which extremize the action. Such extremal worldsheets are often called classical 
trajectories, classical motions, or instantons. 

We thus have: 

(boundaries) "= " ^ g-Action[x{a,T)]_ ^2-24) 

instantons 

It just remains to count how many instantons there are, with given boundary condi- 
tions. 

Typically, we shall require not only that boundaries are fixed, but also that the 
topology of worldsheets be fixed. Concerning boundary conditions, it is certainly pos- 
sible to imagine a very large set of possible boundary conditions (modulo conformal 
reparametrizations again), but we shall consider specific boundary conditions which 
can be classified by some moduli and quantum numbers. In other words, we shall 
consider only certain types of boundary conditions, often called branes, which can be 
parametrized by a finite number of complex variables referred to as open moduli. 

Finally, that defines a function: 

Wt,'\z,,...,zk) (2-25) 

which is the amplitude counting worldsheets of genus g, with k boundaries parametrized 
by k open moduli Zi, . . . ,Zk- 

This paragraph was only a very sketchy and imprecise introduction to string theory. 

Our goal in this article, is to try to understand why those string theory amplitudes 
are the same as those computed by the symplectic invariants and topological recursion 
of [23]. 
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3 Integrability 



Here, we don't assume to have any string theory, instead we are going to construct 
one. Our starting point is a spectral curve as defined by Eq. (12-11) . 

Given a spectral curve, it is always possible to construct a classical integrable 
system (see the reconstruction formula [7]). There is not a unique integrable system 
corresponding to a given spectral curve, but they should all have the same Tau-function. 
Let us choose one of them (this arbitrariness should be linked to the choice of framing 
and background in topological strings, see section 16.51 and in general it is linked to the 
choice of one hamiltonian among the family of commuting hamiltonians) . 

In other words, our starting point is an integrable system, and our goal is to enu- 
merate the classical trajectories given by the equations of motion of this integrable 
system. 

3.1 Action-angle variables and flat coordinates 

Consider a classical integrable system, with a rather arbitrary target space X, with 
coordinates (in a local patch) X = [Xi, . . . , Xo)- Suppose that the motion X{t) is 
a solution of the Hamilton's equations of motion of our integrable theory, or in other 
words it is an extremum of the action (Hamilton- Jacobi equations). 

Then, all classical integrable systems have the property that there is (almost every- 
where), a canonical change of variables {X, dtX) (e*, u) (called action-angle variables) 
which brings the complicated motion X{t) into a linear motion at constant velocity v 
in the Jacobian (see fig. [2]): 

u{t) = u{0) +tv , e = constant. (3-1) 

The Jacobian J is a g-dimensional torus, with some quasi-periodicity properties: 

u = u mod Z9 + rZ0, (3-2) 

where q is the genus of the spectral curve, and r = {7^^} is the Riemann matrix of 
periods of the spectral curve. In fact, it may happen, for some integrable systems, 
that the matrix r be degenerate, so that the periods could become infinite in certain 
directions. In that case, the Jacobian has also non periodic directions, and is a product 
of some power of C times a torus. This situation which may seem non-generic, is 
actually often realized for many examples of interest. Let us ignore it for the moment, 
and assume that degenerate cases can be obtained as limits of the non-degenerate ones. 



13 



Figure 2: Under the action-angle change of coordinates, a comphcated integrable mo- 
tion X{t) in the target space j£, becomes a complex time linear motion u{t) with 
constant velocity in the Jacobian J. In other words, the worldsheet in the target space, 
is a plane (a complex line) in the Jacobian. The Jacobian is a torus, with periodicities, 
and thus the worldsheet can be periodic in some directions. This figure is only an 
"artist view" since the Jacobian should have even real dimension (never dimension 3). 



The complex time evolution of the vector X{t), sweeps a surface embedded into X, 
which we call a worldsheet. In other words, every classical solution of the equations 
of motion, corresponds to a worldsheet. 

The 1-dimensional real curves X{a+iT), a G M at fixed T G M, are called "strings" . 
The worldsheet is indeed the surface swept by a string as T sweeps M. 

After the action-angle change of coordinates, the worldsheet in X is mapped to a 
2-dimensional "plane" (a complex line) in the Jacobian. 

The fiat coordinates on the plane in the Jacobian, can be pulled back to a system 
of fiat coordinates on the worldsheet embedded in X (see fig. [3]). 




Figure 3: The fiat cartesian coordinates on the plane in the Jacobian, provide a system 
of fiat coordinates on the worldsheet. 



Remark 3.1 Another usual formulation of integrability, is the existence of some flat con- 
nection. Here, we see that the flat connection can be realized as the pullback of the parallel 
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transport in the plane, i.e. the pullback of the trivial flat connection in the Jacobian, by the 
action-angle change of coordinates. 

3.2 Boundaries and branes 

There is not a unique choice of fiat coordinates on a plane. Let us see one canonical 
choice, adapted to a choice of specific boundary conditions called Branes. 

In what follows, we wish to enumerate worldsheets having certain topologies, and, 
if we want to have only finite numbers, we need to prescribe some constraints. In 
particular we want to ensure that there is a unique choice of local fiat coordinate t on 
each worldsheet. 

3.2.1 Branes 

For defining the boundary conditions, one first fixes a lattice vector G + rZ^. 
This vector, which we shall call polarization, is a modulus of the boundary which is 
kept fixed from now on. In the following, we consider only boundaries with the same 
polarization. 

Then, we want to consider worldsheets, whose boundaries have the topology of a 
circle and such that in the Jacobian, the boundary is a straight horizontal line parallel 
to v: indeed, the boundary can be a circle only if the straight line in the Jacobian 
is parallel to a periodic direction, i.e. to a period lattice, and we choose it to be the 
polarization v . 

Definition 3.1 D-brane: 

A worldsheet is said to have an D-brane boundary condition with polarization v if 
and only if, in action angle coordinates, the boundary is a straight line in the Jacobian 
parallel (with a real scalar factor) to the polarization v. 

3.2.2 Cannonical choice of fiat coordinate 

Our goal now, is to enumerate worldsheets having D — brane boundary conditions. In 
order to have only a finite number of them, we need to specify some extra conditions. 

We consider worldsheets having D-brane boundaries with a marked point on the 
boundary and with a given "length" parameter / G IR+ (also called "perimeter" of the 
boundary). 

Given a marked point on the boundary and a length parameter I, we choose the 
unique cartesian coordinate t on the plane in the Jacobian, by choosing the origin t = 
at the marked point and the unit such that t = I corresponds to v. 
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The boundary is then the horizontal hne Imt = in the plane, and it is periodic 
of period t = 1. 

Remark 3.2 On a circle of perimeter /, there are / possibilities to mark a point on the 
circle. This means that enumerating worldsheets with a marked point on the boundary, or 
worldsheets without marked points, merely amounts to multiplying by /. 

So, a given marked point and length provide a unique choice of time coordinate on 
the plane, i.e. a unique fiat coordinate on the worldsheet, at least in a vicinity of the 
boundary. 



Figure 4: If we consider worldsheets whose boundary is a circle which is a straight line 
in the Jacobian, and with length /, we choose the unique coordinate t on the plane in 
the Jacobian such that the marked point has coordinate t = and the lattice period 
in the Jacobian has coordinate t = I. That provides a unique canonical choice of flat 
coordinates on the worldsheet, at least in a vicinity of the boundary. 



We call the lines: 

• Imt = constant: horizontal trajectories, 

• Ret = constant: vertical trajectories. 

The boundary is the horizontal trajectory Imt = 0, and the worldsheet is locally 
near the boundary, given by the Poincarre half-plane Imt > 0. 

Remark 3.3 With the same idea, we could also consider "open" worldsheets whose bound- 
aries are vertical trajectories. That would correspond to Von Neumann boundary conditions 
for our branes. We could also enumerate open strings in that framework, but for simplicity, 
we don't do it in this article, and postpone it to a later work. 

4 Decomposition of worldsheets 

We are now interested in "counting" (with Boltzmann weight and symmetry factor) 
all worldsheets, which are orientable connected Riemann surfaces, of a given genus g 




jr 
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(which, once again, has nothing to do with the genus g of the spectral curve), and a 
given number k of boundaries with D-brane boundary conditions of polarization v and 
with a marked point on each boundary of respective lengths li, . . . ,lk- 

Since the image of a worldsheet is the plane of equation 

u{t) = u{0) + jv 

in the Jacobian, counting worldsheets with a given topology, amounts to count all initial 
conditions n(0) G J compatible with the sought boundary conditions and topology. 

So, moduli spaces of worldsheets can be viewed as submanifolds of the Jacobian, 
and are naturally endowed with some measure inherited from the Jacobian. 

Notice, that if the initial condition (the vector 'u(O) in Eq. (13-11) ) is arbitrary, it is 
very likely that the worldsheet will not have the right topology. In fact, it might not 
be compact, it might have infinite genus. Therefore, counting only worldsheets with 
finite genus is very restrictive, we are counting only a very small subset of all possible 
worldsheets. 

This should be related to the fact that we are computing only the perturbative 
expansion of string theory amplitudes, and the non-perturbative part is not captured 
by a genus expansion. We expect that the non-perturbative part introduced in [211 125] 
should take into account those infinite genus worldsheets. 

4.1 Discs 

Suppose that we want to count worldsheets, with brane boundary condition as above, 
having the topology of a disc, i.e. planar g = and only one boundary k = 1. In fact, 
it cannot really be a disc, because the function T = Im t is harmonic on the worldsheet 
and constant on the boundary. That would be impossible on a simply connected 
domain. This means that there must exist at least one singularity of T = Imt inside 
the disc, in other words we have a "punctured" disc. In some sense, this is an infinite 
half-cylinder, but by abuse of language, we shall continue to call it a disc. 

The punctures can sit in some critical submanifolds in the target space, for instance 
the non-compact directions of the target space. Our target space and integrable theory 
may be such that there can exist several kinds of punctures. We shall always need to 
specify which kind of punctured disc we are talking about. 

4.2 Moduli spaces 

Let us consider worldsheets of some genus g, with k brane boundaries of given common 
polarization, and with respective lengths li, . . . ,1^. The topology will be called stable 
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if 

2-2g-k<0. 

Discs {k = l,g = 0) and cylinders {k = 2, g = 0) are not stable. 
We define: 

Definition 4.1 Let M.g,k{PiiW'-,P2ih'i ■ ■ - 'iPki^k) be the set of all oriented connected 
worldsheets of genus g, with k hrane boundaries with marked points and lengths 
h, . . . Jk, sitting on branes labeled pi, . . . ,pk respectively, and quotiented by additions 
of non-singular bare cylinders at the boundaries (see remark below). 

■Mg,k{pi,h]P2j2',---',PkJk) has an orientifold structure, i.e. each worldsheet is 
counted quotiented by its symmetries. 

Let: 

■^t^\piJi;P2,k;- ■ ■■^PkJk) ":=" #-Mg,fc(pi,^i;P2,^2;---;Pfc,/fc), (4-1) 

be the number of elements in Aig^kiPi,h',P2,l'2',---',PkJk) where each worlsheet S is 
counted with a Boltzmann weight q-^^^^^^^C^) coming from the action of our integrable 
system, and with a symmetry factor l/#Aut(S) if it has non trivial automorphisms. 
More precisely: 

„— action (S) 

^ft\plJl■,P2J2■,...■,p,,h) := sAntm - ^^'^^ 

Notice that stable Riemann surfaces (i.e. 2 — 2g — k < 0) have a finite number of 
automorphisms . 

The Boltzmann weight e~^'^*'°'^ depends on the integrable system, i.e. on the moduli 
of the target space, as well as some additional parameters and the moduli of the 
brane boundaries. We will not have to know the weight, and we shall prove that 
the topological recursion holds without having to specify the weight. This weight is 
encoded in the spectral curve or in the disc amplitude. 

Let us explain what we mean by counting worldsheets, "modulo addition of non- 
singular cylinders at the boundary" (see fig ED- This means that if a worldsheet is 
obtained from another, by analytically extending the fiat coordinates near the bound- 
ary, both worldsheets are in the same equivalence class and should be counted only 
once. This means that Aig^k{pi, h]P2, h', ■ ■ ■ ',Pk, h) is locally independent of pi, . ■ . ,Pk- 
It is invariant under small changes of pi, . . . ,pk. But it can change when one of the 
Pi's approaches a special brane, a branching or a puncture. 

The purpose of the next sections is to show how to compute those "weighted num- 
bers of worldsheets" in terms of the spectral curve of the integrable system. We will 
see that they can be computed recursively by the topological recursion. 
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Figure 5: Worldsheets ending on brane pi or p[ are considered equivalent, and are 
counted only once in A4i^2- But they are not equivalent to worldsheets ending on ai. 



4.3 Branchings 

Let us consider a worldsheet in J^g^k{pi, h]P2, h] ■ ■ ■ ]Pk, h), with 2 — 2g — k < and 
A; > 0. 

Notice that if 2 — 2g — k < 0, the worldsheet is neither a cylinder nor a disc. It is 
locally a half-cylinder near its boundaries, but it cannot keep the topology of a cylinder 
under time evolution, i.e. it cannot be globally mapped to a plane bijectively. There 
must be some time at which the flat coordinate becomes singular. 

Consider one of the boundaries (let us say the first one pi,li) with its marked point 
and length li, and choose the unique flat coordinate t defined as above in the vicinity 
of that boundary. 

For small times Imt, the worldsheet is a cylinder, and horizontal trajectories are 
circles winding around the worldsheet. However, since the worldsheet is not globally 
a cylinder, there must exist a time at which the horizontal trajectory is no longer a 
circle. 

Several situations can occur: 

• The horizontal trajectories can hit another boundary. This would lead to a 
worldsheet with the topology of a cylinder since all boundaries are parallel in the 
Jacobian, that is to say g = and k = 2 which is not considered here. 

• The horizontal trajectories can hit a puncture. This would mean that our world- 
sheet would be a disc, i.e. g = and k = 1, which is not either the situation we 
wish to consider here. 

• They can hit a singularity where the horizontal trajectory gets pinched and splits 
into two (or more) connected components, the generic situation corresponding to 
two components (see fig. [6]). The circle at time Imt < tc, splits into a "figure of 
8" at tc- After tc, the worldsheet splits into two half- cylinders. 
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• In fact, just by going backwards in time, we see that there are also singularities 
at which the horizontal trajectory becomes a half of a figure of 8, where another 
half-cylinder could join. Again, at t = tc, the horizontal trajectory is made of 
two circles, i.e. a "figure of 8". 

Let us call "branch points" Oj, the points in the moduli space of branes (whatever it 
is) at which such branchings may occur. From now on, we assume that our integrable 
system be such that we have only a finite number of 





Figure 6: The horizontal trajectory hits a singularity for the first time at time t^, at 
which the fiat coordinate is ill-defined and the circle splits into a figure of " 8" . Two 
possibilities may occur: after the cylinder splits into two half-cylinders, or it merges 
with another cylinder to make a bigger cylinder. Remark that the second possibility is 
the same as the first one, under time reversal. In both cases, the horizontal trajectory 
at time Imt^ is a figure of "8". Situations where more than two cylinders join, or split, 
are non generic, and should bring a vanishing contribution to the generating function 
counting worldsheets. 



Definition 4.2 We shall call a "bare propagator" a piece of an open worldsheet, < 
Imt < Imtc which is a cylinder, where the flat coordinate is globally defined, and all 
horizontal trajectories are circles, and such that at t = t^, the horizontal trajectory is 
a figure of 8 (see fig. [^. 

We denote such a bare propagator: 

S{poJo;aiJ) (4-3) 

where the horizontal trajectory Imt = is on brane po, with length Iq, and the horizontal 
trajectory Imt = Imtc is on the critical brane ai, with length I = lo- 

The main property, is that beyond time tc, the horizontal trajectories should be 
conformally continued into two half cylinders. 

^Remark that the branch points depend on the choice of integrable system and polarization. 
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Figure 7: The bare propagator: its interior is a cylinder, bounded by horizontal trajec- 
tories. One side is a circle, the other side is made of two circles glued at a branchpoint. 
Either the circle of length Iq = I is pinched into two circles I = I' + {I — I') , or it is one 
half of a pinched circle of total size I" = I + V . 



4.4 Recursive decomposition of the worldsheet into propaga- 
tors and cylinders 

Consider (7 > and A; > such that 2 — 2(7 — (/c + 1) < 0, and consider a worldsheet 
S G A^g,fc+i(po5 ^0; P17 h]'P2^ h', ■ ■ ■ Pk, h)- Consider its first boundary, of length Iq, with 
a marked point, ending on brane po, and consider the unique flat coordinate t on the 
worldsheet near this boundary. 

Since the worldsheet doesn't have globally the topology of a disc or cylinder, there 
must exist a smallest time t^, Imtc > 0, at which the flat coordinate becomes ill defined 
and a branching occurs. In other words, there exists a, generically unique, time tc at 
which we reach some branchpoint Oj, and therefore the worldsheet S contains a bare 
propagator S{po, /q; Oj, /). We emphasize that this propagator is uniquely defined. 

Let us call S' the worldsheet obtained by removing the bare propagator from E: 

S' = S\S(po,/o;a„0. (4-4) 

S' has again brane boundary conditions (boundaries are indeed horizontal trajec- 
tories parallel to v), it has k + 2 boundaries since one of the boundaries is split, and 
S' is either connected or disconnected: 

• If S' is connected, it is clear that it belongs to either Aig^i^k+2{ci-i,^'] — 
Pi, /i; P2, ^2; • • • ; Pfc, k) or A<g„i,fc+2(ai, Oi, / + pi, /i; pa, h\---\ Pk, h), i-e. in both 

cases to 7Vlg-i,fe+2(a.j, \l'\;ai, {l - l'\;Pi, h;P2, k] ■ ■ ■ ]Pk, h) for some /' G M. 

• If S' = U S'^ is disconnected, the two connected parts belong to 
Aih,i+#i{cii, I) and A4h',i+#i'{0'i, \^ ~ ^'\', I') for some h, h', /, /' such that h + h' = g 
and / l+l /' = {pi, h; . . . ; Pk, h}, and /' G M. 

When it is disconnected, it may happen that one of the two connected components, 
let us say is a punctured disc G A4o^i{ai, |/ — /'|), and the other connected component 
then belongs to A4g^k+i{cLiA^'\iPijk',P2j2',---',PkJk), like S itself (in particular, 
the other connected component can't be a disc). In that case we may redo the same 
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Figure 8: Consider a worldsheet E G M.g^k+i{po, Pi, h] ■ ■ ■ ] Pk, h) which is not a disc 
or a cyhnder. Start from the boundary {po,lo), consider the horizontal trajectories 
defined from that boundary. There must exist a smallest time tc at which the hori- 
zontal trajectory stops being a circle. The piece of surface < Imt < Imt^ is a bare 
propagator. If we remove the bare propagator and the figure of 8 critical trajectory 
from S, we get a worldsheet S'. S' is either connected or disconnected. 



thing on S'_: start from the boundary, until we reach a branchpoint, and remove the 
corresponding bare propagator. We can do that recursively, until none of the connected 
components is a disc. 

It is thus more convenient to define a " renormalized propagator" which may include 
an arbitrary number of discs glued. It is defined by the property (see fig. Mj 

S{poJo;ai,l) 
= d{lo - I) S{po,lo;ai,l) 

/oo ^ 
dl' S{po, Iq; aj, Iq) U S{aj, |/'|; Oj, /) U Disc(aj, |/o — /'|). 

(4-5) 
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Figure 9: The renormalized propagator, is obtained by following the horizontal trajec- 
tories from the first boundary. Each time a critical trajectory is met, the surface may 
split into two disconnected parts. We recursively do that until none of the connected 
components is a disc. The renormalized propagator thus contains a certain number of 
bare propagators and discs. It ends at a critical trajectory. 



Removing the renormalized propagator from a worldsheet S, means removing a 
propagator from the first boundary, and if one of the connected components is a disc, 
then remove again a propagator from the other connected component until none of the 
connected components is a disc: 

S" = S\S(po,/o;a„/). (4-6) 

Notice that, since none of the connected components of S" is a disc, then, each con- 
nected component of S" has a Euler characteristics strictly larger than that of S, and 
therefore after repeating this procedure a finite number of times, we arrive to only 
propagators and cylinders. This is a topological recursion. 

Therefore we can decompose any worldsheet S, in a generically unique way, into a 
finite number of renormalized propagators, and cylinders. 

Finally we have the following (orientifold) bijection between moduli spaces: 

Mg,k+i{poJo;piJi;P2j2; ■ ■ ■ ;PkJk) 

POO POO _ 

- ^ / dl dl' Ms{poJo;aiJ) X Mg^i^k+2{ai,\l'\;ai,\l - I'\;puh;p2,l2; ■ ■ ■;PkJk) 

^'h+h'=g;i^i'={pi,h;-;Pk,ik} Mh,i+#i{ai, |/'|; /) x Mh',i+#i'{ai, \l - (4-7) 

where U' means that we exclude discs i.e. (/i, /) = (0,0) and {h',I') = (0,0), and 
■Ms{po, Iq] O'i, is the moduli space of all renormalized propagators with one boundary 
of length Iq on the brane Pq, and the other boundary on brane a^, of length /. 

This recursive decomposition of moduli spaces of worldsheets is very similar to the 
recursive structure of the topological recursion Eq. (12-201) . In the next sections, we 
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shall show that the generating functions of the volume of these moduli spaces indeed 
satisfy the topological recursion. 

4.5 Skeleton graph of the worldsheet 

Another consequence of that decomposition, is that we can associate a graph to any 
worldsheet. 

Indeed, chose the first boundary and remove renormalized propagators, i.e. remove 
propagators until none of the two connected components is a disc. Then, draw the 
splitting horizontal trajectory, that is the boundaries of the renormalized propagators, 
on the worldsheet, and proceed recursively, until it remains only cylinders. 

We have thus drawn some dividing circles on each world sheet. On each renor- 
malized propagator, let us draw an arrowed line from the marked point on the circle 
boundary, to the branchpoint. On each cylinder, let us draw an unoriented line between 
the 2 marked points. 



We obtain a graph drawn on each worldsheet, whose vertices are labeled by branch- 
points. These are exactly the graphs of [23]. Those graphs have 2g — 2 + k vertices, 
2g — 2 + k arrowed edges, k — 1 external non-arrowed external edges, g internal non- 
arrowed edges, forming g loops, and such that the 2g — 2 + k arrowed lines form a 
tree rooted at the first boundary and going through all vertices. The edging is also 
constrained by the following rule: non-arrowed lines can only connect two vertices if 
one is the descendent of the other along the arrowed tree, see [20l [23] . 

Remark 4.1 The graph obtained depends on a choice of a "first boundary", another choice 
could lead to another graph. 

One can see that the dual of the graph (see fig. [TO!) . obtained by drawing circles 
dividing every cylinder into two half-cylinders, provides a pant decomposition of the 
worldsheet. But contrarily to the usual Teichmiiller spaces approach, here, thanks to 
our integrable system, we have for each worldsheet and choice of "first boundary", a 
unique pant decomposition. We don't have to consider a quotient by the mapping class 
group. In some sense we have already chosen one canonical representant in each class. 

Remark 4.2 Different worldsheets in Mg^k can have different graphs. This allows to de- 
compose the moduli space A^g,fc into a finite number of cells labeled by graphs. Somehow 
this is in the same spirit as the Strebel-foliation used by Kontsevich, or the Teichmiiller pant 
decomposition. 
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Figure 10: Example: the graph on a worldsheet of M.i^2- Each worldsheet has a unique 
graph, once we have chosen an entrance boundary. Notice that the dual (blue circles in 
the middle of each cylinder), gives a canonical pant decomposition of the worldsheet. 




Figure 11: The two possible graphs obtained for worldsheet s in Aii^2- 



5 Topological recursion for the amplitudes 

In this section, using standard methods of combinatorics, and in particular Laplace 
transforms, we translate the recursive decomposition of our moduli spaces into relations 
for the generating functions N'lf\pi, h',P2j2',---',Pk-,h)- 

In the preceding section, we showed that integrability allows to get a unique foliation 
of every worldsheet through the use of a flat connection built from the integrability 
of the considered system. This allows to get a cell decomposition of the moduli space 
of surfaces labeled by the graphs of [201 123] described in section 14. 5[ Moreover, this 
gives a bijective procedure to build every worldsheet by gluing discs and cylinders. In 
this section, we translate this bijection into a recursive relation among the amplitudes 
■^t\piJi]P2,h', ■ ■ ■ ',PkJk) which are the generating functions counting the elements 
of Mg^k (with Boltzmann weights and symmetry factors). 
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5.1 Disc amplitude 

For the moment, we shall not explain how to compute the disc amplitudes J\fi^\pi, li). 
We shall assume those numbers to be given. The way they depend on the integrable 
system, or how they are related among themselves, will be explained later in section 
16.41 Let us just mention that they are "difficult" to compute. 

We also emphasize that what we call a disc is really a " renormalized disc", i.e. it 
may contain several punctures and branchings, the flat coordinate needs not be globally 
well defined on it. Renormalized discs can be obtained by gluing propagators and bare 
discs recursively, in the same way we defined the renormalized propagator. But since 
we don't know yet the generating function of bare discs, we shall not perform that 
construction here. 

For an brane, we define the generating function as a function of a formal complex 
variable z G C, by a Laplace transform: 

/•oo 

Wi°\ai,z)= dl e-'' Uf'\ai,l). (5-1) 

Notice that M^\ai^l) counts discs with a marked point, and therefore: 

ATf (a„0 = /A/-r(a„/) (5-2) 

where J\f^^\ai,l) counts discs without marked points. This implies that wl^\ai^ z) is 
a derivative with respect to z, and this is why Wi^\ai, z) dz should be thought of as a 
differential form. 

From now on, we shall assume that our integrable system is "regular", i.e. such 
that wf^\ai^ z) is analytical near 2; = 0, and that W^\ai^ z) + w'f\ai^ —z) has only a 
double zero, and not a higher order zero. We say that our integrable system is critical 
when we have higher order zeroes. 

5.2 Propagator amplitude 

The bare propagator means the non-renormalized propagator, i.e. a piece of a world- 
sheet, where the fiat coordinate is globally well defined. Let us denote its amplitude 
by: K{j)Q, Iq; ai, I). Because there is no singularity, the lengths are conserved, and thus 
the amplitude can be non zero only if /q = /. 

The renormalized propagator is obtained by gluing discs to one of the 2 boundaries, 
i.e. its amplitude K{pQ,lo;ai,l) satisfies the relation (see figE]): 
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S{1 - lo) K{po,lo;ai,l) 

/oo 
dl' K{poJo;a,Jo) ATf (a„ |/o -/'|) K{a,, |/'|;a„/). 
■oo 



(5-3) 



This relation doesn't determine K or K, and we shall see later how to determine 
them. However, it allows to express K in terms of K and 

In Laplace transform we define the generating function: 

POO 

Kpo,io{auz)= dl e'^^ Kp^^i^.aui- (5-4) 

^0 

5.3 Annulus and propagator amplitude 

We call annulus or cylinder the elements of A^o,2- 

Let us consider a bare cylinder, i.e. a worldsheet starting on a brane pi with length 
li, and ending on p2 with length 12- "Bare" means that the flat coordinate is globally 
well defined on the cylinder and doesn't encounter any singularity between the two 
boundaries. 

It is obvious then, that all horizontal trajectories have the same length and thus 
h = I2, i.e. 

77fipuh;p2,l2)=0 ifhy^k 

for any pi,P2- Also, since we count worldsheet s modulo gluing a bare cylinder, we see 
that if the boundaries pi and p2 are close enough, all cylinders are equivalent. 

If pi and P2 are close enough, it is clear that there is only one possible cylinder, and 
there is / possibilities of choosing a marked point on the second boundary, so that: 

A/'2°^(pi,/;P2,0 = ^• 

Then, consider the renormalized cylinders. This means that between branes pi and 
P2, the flat coordinate t may encounter branchings where one branch ends on a disc. 
In other words we may glue many discs recursively. The process of gluing discs at 
branching is already included in the renormalized propagator (see fig [T2l) . so that we 
have: 

ATf (pi, h;p2, h) = 5{h - I2) ATf (pi, /i; p2, k) 

poo 

+2J2 dl' K{p^, k; aj, I' - k) Af}^\a„ I') Jffia,, l2]P2. k] 
+2 V dl' K{p,,k;a„k-n M'\a„l') Jff{a^MP2M 
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poo 

+2 V / dl' K{p,, h; a„ I' + h) U'f\a,, I') AT^ (a„ W,p2, h) 
6{li - /2) Ar2°\pi, /i;p2, k) 

/oo 
dl' K{p,, h; a,, \k - I'D AT^ (a,-, k;p2, k) 

-oo 



(5-5) 



If the second brane p2 = ai is at a branchpoint, this means the cyhnder ends exactly 
where the flat coordinate degenerates. The cylinder is thus exactly the same as the 
propagator, except that we need to glue a disc at the other half-boundary, so that the 
boundary has the topology of a circle. That translates to: 



^^f2°\pl,k■,ai,k 



dl' K{puluai,\k-l'\)M'\a,,\l'\) (5-6) 



Notice that we have a l/k in the left hand side, and a 1/2 (in fact 2 times 1/2) in the 
right hand side, because of symmetry factors of marking a point on the boundary. 




Notice that this last relationship, together with Eq. (15-51) implies: 

N'!f\pi,k]P2,l2) = S{k - k)J^2\piJi;P2Ji) 

+ M!f\pi,k]aj,k) Jf2\(^j^k]V2^2)- 

which just means that gluing two cylinders gives a cylindei0. This is the self- 
reproducing property of cylinders. 

We define the generating functions: 

POO 

B{po, no; ai, z) = dl e'""^ M!f\pQ, k; ai, I) (5-7) 



Eq. (15-61) translates into 

B{pQ,k;ai,z') dz' 



^Since we consider marked points on the boundary, we need to divide by I2 in order to forget the 
marking on the intermediate boundary. 
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~ dlo 



Jo 

oo /"OO 



dl / dl' e-'^^-^'^ K{j>Q,lo]ai,l)M['^\ai,l') 

dl / d/'e+"(''-') K(po,^o;ai,0-^r(«i,0 

+ / dl I dl'e-'^^'^^'^ K{po,lo■,ai,l)^^}°\a^J') 

dl I dl'e-'\e'^' +e-'^') K{pQ,lo;ai,l)U[^\ai,l') 



Ji 

oo /"OO 











(0), 



I.e. 



K{po,lo;ai,z) 



(5- 



(5-9) 



{Wi'\ai,z) + w;'\a,,-z)) 
This relationship is merely the combinatoric relation illustrated by fig [T2|, and is to be 
compared with Eq. (12-81) . 

5.4 Topological recursion 

The bijective procedure of section 14.41 tells us that moduli spaces of stable topologies 
can be decomposed recursively. The bijection Eq. ( 14- 7p clearly translates into the 
following relation among amplitudes: 

■^Si{poJo;piJi; ■ ■■■,Pk,lk) 

poo poo i- 

= X] / '^W dl' K{pQ,lQ-,ai,\l + l'\) Afl:l~^\ai,l-,ai,\l'\-,pi,li,p2,l2-, ■ ■ -^PkJk) 



+ 



A/'S',,(''i.'i-')-Vi+'«'(''..l''li-'') 



h+/i'=g;/l±l/' = {pi,Zi;...;pfc,«fc} 

OO poo poo 

dz I dl" I dl / dl' e-^('"+'+'') 







oo J — OO 



(5-10) 



K{po,lo;ai,l") J^k^+2^\ai, \l\;ai, \l'\;pi,li;p2,l2; ■ ■ ■■,Pk,lk) 

+ E M[%M^^l\■J)KlI'Ml'\■J') 

h+h'=g;HHl' = {pi,h;...;pk,lk} 



where Yl' means that we exclude {h, I) = (0, 0) and {h', I') = (0, 0). 
With the Laplace transforms generating functions 

poo poo 

Wlf\ai^,zi; . . .■,ai^,Zk) = / ... / dli...dlk Afl:^\ai^,li; . . . ; ai^Jk) Y\. ^~ 

Jo Jo j^-^ 



(5-11) 
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the recursion relation can be rewritten 

jig) 

k+l 



W^i+i(aio' ^o; Oj^, 2:1; . . . ; ai^^Zk) 



= y^^TT— dzK{aio,zo]ai,z) wlf,2\ai,z]ai,-z]ai^,zi]...]ai^^,Zk) 

+ W?i%,(a„^;/)iy/fi,,(a„-z;/')]. (5-12) 

/i+/i'=g;/l+)/'={ai^,2i;...;ai^,2;fc} 

By an easy recursion, one sees that the integration contour over z can be deformed 
into a circle surrounding the pole at z = 0, i.e. 

ERes dz K{aig,Zo;ai,z) Wlf_r2^\ai, z; ai, -z; ai^, Zi, . . . ]ai Zk) 
2^0 L 



+ W[%j{a,,z-l)w[%,{a,,-z-l')\. (5-13) 

ft.+ft.'=g;/ai''={ai-^ ,zr,...]ai^,Zk} 

This is the topological recursion written in terms of a local coordinate z near each 
branchpoint. 

The only thing we need to do, in order to fully recover the topological recursion of 
[23], is rewrite all those residues formula in terms of intrinsic variables on the spectral 
curve, rather than local coordinates. 

5.5 Closing a boundary 

Consider 2 — 2g — n < Q and a worldsheet in M.g^n{pi^ h] ■ ■ ■', Pn, In)- As we have seen, 
it is obtained by gluing along critical horizontal trajectories 2g — 2 + n renormalized 
propagators, and g + n — 1 renormalized cylinders. 

Remember that a renormalized propagator S'(pi, /i; Oj, /) is a cylinder whose second 
boundary ends on a critical horizontal trajectory. The first boundary ends on brane 
Pi- 

However, in the topological recursion, we need to consider also renormalized prop- 
agators, whose both ends are on critical trajectories S{aj, a^, /). In order for the first 
boundary to have the topology of a circle so that we can glue it to another propagator, 
we need to glue a disc on the second connected component of the first boundary. 

This means that every internal renormalized propagator must contain a disc. Notice 
that for the cylinders, we have the possibility of having the bare cylinder, which contains 
no disc. This means that generically, a worldsheet of M.g,n contains 2g + n — ?> discs. 

If the boundaries pi,...,pn are themselves critical branes aj^,...,aj^, the first 
propagator starting on a^^ must also contain a disc, and therefore, a worldsheet of 
■M.g,n{o,h, h', ■ ■ ■ o,i„, In) coutaius 2g + fi — 2 discs. 
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Figure 13: Example of a worldsheet in 7^2,1- It can be decomposed into 3 propagators 
and 2 cylinders. Since the initial boundary of each propagator is a critical trajectory, we 
need to close the second component of the "8" by a disc. In other words the worldsheet 
must contain 3 discs. By cutting out one disc, we get a worldsheet in A^2,2- 



Consider a worldsheet S G M.g^n{(iii, h] ■ ■ ■ ] o,i„, In) (see fig fT3|) . Choose one among 
the 2(7 + 77,-2 critical trajectories at the initial end of an internal propagator, and cut 
the worldsheet along that trajectory. It gives 

S = S' U Disc(aj, /). , S' G Mg.a+i{ai^,li] ai^, /„; a^, /). (5-14) 

Since there are 2g + n — 2 ways of doing that, we have a.2g + n — 2^1 application 
from M.g^n+1 ■^g,n- This implies: 

{2g+n-2)Ulf\ai^,li]...]ai^,ln) = / y TV^f ^(ai,, /i; . . . ; a^^, /„; a„ /) Mi^\ai,l) 

i "^0 

(5-15) 

(we need to divide by / because otherwise the marked point is marked twice, once in 
J^n+i and once in 

In Laplace transform that gives: 

{2g+n-2) w!;{'\ai^, Zi] ai„, z^) = Res wjf^^{ai,, Zi, ai„, z^, ai, z) $(0^, z) dz 

(5-16) 

where 

rf$(ai, z)/dz = Wi^\ai, z). (5-17) 
This relationship was also derived in [23] as a consequence of the topological recursion. 

5.6 Closed surfaces 

Worldsheets belonging to M.gfl have no boundary. However, we shall assume that 
they are defined also with respect to the same polarization v. This means, that any 
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worldsheet S G Mg.o, is also a plane parallel to v in the Jacobian. We may choose any 
point on E, and draw the horizontal trajectory going through it. 

So, let us choose a point on S, and let us cut the worldsheet along the horizontal 
trajectory through that point. The resulting worldsheet S' maybe disconnected or not. 
It has two boundaries with brane boundary conditions, and thus S' belongs either to 

Mg-1,2 or Mg',1 X Mg-g',1. (5-18) 

The results of previous section imply that S' can be decomposed into propagators and 
cylinders, and must have 2g — 2 discs (we assume g >2). 

Let us choose one of the discs of S, and choose a point on the boundary of that 
disc. Let us then redefine S' as S cut along the horizontal trajectory going through 
that point. 

In that case we have 

eMo,i{aiJ) X Mg,i{a,,l). (5-19) 

In other words, every worldsheet S G Mg,o can be obtained by gluing a disc to the 
boundary of a worldsheet G -Mg,!, which can be decomposed into 2g — 1 propagators 
and g cylinders. 

Therefore every worldsheet S G Mg,o can be decomposed into a disc, 2g — 1 prop- 
agators, and g cylinders. 

This decomposition is not unique, it can be done for any of the 2g — 2 discs, this 
means that the decomposition is not bijective, but is — 2 — > 1. 
We thus obtain: 

W^o'^ = Fg = V Res Wt\ai, z) $(a„ z) dz. (5-20) 

^ i 

This is precisely how F^'s are defined in |23j . 

6 Reconstruction of the spectral curve 

So far, we have defined generating functions w'^'^ as formal series of complex formal 
variables Zi-s,. Here, we shall glue all patches in order to get generating functions 
globally defined on a Riemann surface C. Then, afterwards, we shall show that this 
Riemann surface C, is actually the same as the underlying Riemann surface of our 
starting point spectral curve C = C. 
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Figure 14: Example of a worldsheet in A^2,o- Chose an arbitrary point on the world- 
sheet, and cut the worldsheet along the horizontal trajectory going through that point. 
One may get a worldsheet in Aii^2 oi Aii^i x Aii^i. Assume that we are in the A4i^2 
case. A worldsheet in A^i,2, is obtained by gluing two propagators and 2 cylinders. The 
internal propagator has its initial boundary on a critical trajectory, and its boundary 
must be a circle, so that we need to close one half of it by a disk. This means that there 
is a disc on our initial worldsheet. Now, cut the initial worldsheet along the critical 
trajectory of the disc boundary. One gets a decomposition of our initial worldsheet 
into a disc and a worldsheet in Ai2,i- 



6.1 The disc amplitude 

Consider a Hurwitz space {C,x), given by a Riemann surface C of genus g (the same 
genus as C defining the spectral curve of the integrable system, i.e. the dimension of 
the Jacobian), and a projection x : C ^ CP^, with as many simple ramification points 
as the Cj's. 

Assume here that dx is a meromorphic form, whose zeroes are simple, and are 
labeled by the a^'s. 

The degree of dx (the number of poles with multiplicities) is then #{aj} — (2g — 2). 

Notice that if C e C is a point on C near Oj, then z = \/ x(C) — x[ai) is a local 
coordinate on C near Oj, and we have: 

x(C) = x(a,) + z^. (6-1) 

The two branches C and C coming together at the ramification point, correspond to z 
and —z respectively , i.e. 

x{Q=x{Q , ^Jx{l)-x{a>) = -^x{(:)-x{a,). (6-2) 
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Now, let us try to find an analytical function defined on an open domain of C, 
containing all branchpoints, and such that in the vicinity of branchpoint Cj, we have 
in the local coordinate z = \/ x — x(aj): 

2^ m) - m) = Vx(C) - x(a,)) + (a„ - ^£(0 - x(a,)) . (6-3) 



The existence of such a function globally defined of C is not obvious, so, let us 
assume that we are in a situation where it does existl^. 

Then, it is clear that if such a function does exist, it is not unique. Indeed, one 
may add to y any rational function of x: locally near any branchpoints we can add 
any even function oi z = \/ x ~ x{ai). In other words, y is not unique, it is defined up 
to an additive rational function of x. 

So, let us assume that we have chosen one function y on C. 

6.2 The 2-point function 

The cylinder amplitude should satisfy the self-reproducing property, i.e. the fact that 
gluing a cylinder to a cylinder, is again a cylinder. This should hold only if gluings are 
performed away from singularities (branchings or punctures). 
In other words, for any pi,p2,P3, and any li and l2'- 

/CO 7J 
J ATf (pi,/i;p3,0A/'f (P3,/;P2,/2) (6-4) 

(we divide by / so that the marked point is not counted twice). 

If we want to translate that into a global property for the generating function, we 
write: 

W2°\pi,zi;p2,Z2) = Res W2°\pi, Zi;p2, z) dz / W2^\pi, z';p2, Z2) dz' . (6-5) 

z^zi J 

which implies that near zi = Z2 we must have: 

W^°\pi, Zi]P2, Z2) ~ 7 ^ — (6-6) 

[zi - Z2Y 

in any local variable 2, and for any pi,p2 close enough so that we can choose the same 
local variable for zi and Z2. 



®As it is pointed out in the next section, in physics, one goes the other way round: from a problem 
in physics or mathematics, one derives such an integrable system by the computation of the simplest 
observables whose generating function satisfy an equation defining the spectral curve. The existence 
of such a globally defined function is thus ensured from the beginning. 
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On a Riemann surface C, there exists a Bergman kernel, i.e. a symmetric 2-form 
B{Ci,C2), having a double pole at = (2, with vanishing residue, and no other pole, 
and normalized so that: 

^1-22 (2;(Ci) - z{C,2)Y 

where stands for a point of C, and z{C^i) G C is any local coordinate on C. 

Such a form is the definition of the "Bergman kernel", see Eq. (12-21) . One ambiguity 
remains. In order to be uniquely defined, the Bergman kernel is defined normalized 
on a chosen symplectic homology basis of non-contractible cycles on C, see Eq. (12-41) . 
What we can say, is that given a curve C, and functions x and y, there is not in general 
a unique choice of a Bergman kernel. We need another data, independent from the 
disc amplitude, and this data is a g x g symmetric complex matrix k, as in Eq. (12-41) . 

Let us assume that we can choose a Bergman kernel, globally defined on C, such 



that near every two branchpoints Oj, we have (where Zii^Q = yx{(i) — x{ai) ): 

W^f (a„ ^i(Ci); a„ ^,((2)) rf^i(Ci) dz,{C2) = 5(Ci, C2). (6-8) 

6.3 Higher topology amplitudes 

Then, the topological recursion implies by an easy recursion, that all generating func- 
tions of the type 

W^^\ai^,zi; 0^2, 22; • • • ; ai^,Zk) (6-9) 
can be defined globally as meromorphic differential forms on the curve C: 

wj;^^\ai^,Zi,{Ci);ai^,Zi^{C2); ■ ■ . ; ai^, Zi^{Ck)) dzi^{Ci) ■ ■■dzi^{Ck) 
= wt\C^,...,Ck). (6-10) 



Since at each step in the recursion, residues are computed at branchpoints, this implies 
that the forms PF^^^'s have poles only at branchpoints. 

Therefore, we have obtained that the generating functions "counting" worldsheets 
in A4^k \ are the correlators of [23] obtained by the topological recursion from the 
spectral curve (C,x, y). 

Similarly, generating functions Fg of worldsheets with no boundary, are the sym- 
plectic invariants Fg of the spectral curve {C,x,y). 
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6.4 Reconstructing the integrable system 

It was claimed in [23], and made precise in [25] that, out of the symplectic invariants 
and correlators of a spectral curve (C,x, y), it is possible to construct a formal Tau- 
function (formal function of Qs), satisfying Hirota's equations: 

Tig,) = exp {Y.9s'-' F,{C, x, y)^ (l + g^iO' F[ + ^ 6'" F^") + . . . ) (6-11) 

(see the exact expression in [25]). 

This Tau function defines an integrable system, whose classical limit at g, 0, is 
a classical integrable system of spectral curve {C,x,y). 

Therefore, starting from some integrable system, whose classical limit has a spectral 
curve {C,x,y), we have, by using the flat connection, defined a "topological string 
theory" , whose amplitudes define themselves another integrable system whose classical 
spectral curve is {C,x,y). Those two integrable systems encode the same information, 
they are dual to one another, and therefore they should have the same spectral curve 
up to symplectic transformations. 

In other words, up to symplectic transformations, the spectral curve (C, x, y) is the 
spectral curve of the initial integrable system. In particular the underlying Riemann 
surface is the same: 

(C, f , y) = (C, X, y) modulo symplectic transformations. (6-12) 

Therefore, a posteriori, we determine the disc amplitudes by W^'^ = ydx, and the 
cylinder amplitudes as the Bergman kernel on C. 

6.5 Ambiguity of the construction 

Let us notice that we have made some arbitrary choices at some points. 
Framing and choice of an integrable system 

First, given a spectral curve, we have chosen a realization of a classical integrable 
system attached to it. This choice is not unique. In particular, given a curve C, we 
have chosen a projection x : C —>■ CP^, to obtain a Hurwitz space, and we have chosen a 
function y (for instance we have seen that we have the freedom to add to y any rational 
function of x). The Tau-function (and the -Fg's) are unchanged if we choose another 
representant of the same integrable system, in particular we may change [x, y) (x, y) 
such that dx A dy = dx A dy, without changing the F^'s and the Tau-function. But 
when doing that, we change the open amplitudes Wn^^s. 
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This ambiguity can be thought of as a choice of "framing". 

For example, in the context of topological B strings, the mirror spectral curve is of 
the form 



where is a polynomial. For any integer /, changing y ^ y + f x doesn't change the 
Fg's and wlf^^s, and it corresponds to 



which is a well known framing transformation in topological B strings. 
Polarization and modularity 

Then, we have chosen a polarization vector v in the lattice + rZ^. This polarization 
vector can be viewed as a characteristics in the Jacobian. It can be linked to a choice 
of a symplectic basis of non-contractible cycles Ai fl Bj = 6ij on the curve C. Indeed, 
a modular transformations (i.e. a change of symplectic basis {Ai, Sj)), is equivalent to 
a SlgCZ) transformation of the lattice, and can change v to any other lattice vector. 
Notice that we have the same ambiguity in the choice of a Bergman kernel on C. In 
other words, the choice of the Bergman kernel, should be linked to the polarization 
vector V. 

The open amplitudes Wn^^s and closed amplitudes -Fg's, depend explicitly on the 
choice of polarization, i,e, they are not invariant under modular transformations. It 
was proved in [231 122], that modular transformations of the Fg's and Wn^'s, obey the 
formalism of BCOV [8], and are given by the diagrammatic rules of [2]. 

BCOV |8], and [2] noticed that those modular changes of the F^'s can be canceled 
by adding some non-holomorphic terms to them. It was long debated what the role 
of those non-holomorphic terms could be, and in particular, there is no such non- 
holomorphic terms in the Chern-Simons field theory which is supposed to be dual to 
the topological string B-model. 

Recently it was discovered [25] that the modular changes of the F^'s can also be can- 
celed by some holomorphic non-perturbative terms. In other words, the holomorphic 
anomaly is merely an artifact of perturbative expansion. 

In fact those non-perturbative terms are essential to make the whole partition func- 
tion satisfy Hirota equations and be the Tau-function of an integrable system, and we 
used them in section 16.41 above. 



i7(e",e^) = 



(6-13) 



ev ^ (e") 



(6-14) 
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Non-perturbative part 

The whole string theory partition function thus contains a perturbative part, given by 
the Fg^s, which count open strings of finite genus g with a given polarization, and a 
non-perturbative part, which restores modular invariance, background independence 
and integrability. It would be interesting to understand what worldsheets the non- 
perturbative terms count. A guess is that they count worldsheets which are non- 
compact Riemann surfaces, of infinite genus. 



7 Other expansions, Lagrange inversion 

In the previous section, we have constructed some string amplitudes w!f^ which are 
meromorphic forms intrinsically defined on the curve C. They were constructed from 
a spectral curve S = {C,x,y). 

By construction, when we expand w!f\(i, . . . , („) as Laurent series in terms of the 
local variables Zi = \/x{(i) — x{a) near a branchpoint a, the coefficients of the expan- 
sion of w!i'^ are Laplace transforms of the generating functions counting worldsheets 
of topology {g,n) and with boundaries of given lengths: 

where 

h\ , , ..fhji 

Notice that, when 2 — 2g — n < 0, from the general property of the topological recursion 
we have that A^^],^. = ii ki > Qg — 6 + 2n, but the fcj's can take negative values 
down to — oo. Notice also, that if /cj > then ki must be even, and for ki < 0, there is 
no parity restriction. 

For 2 — 2g — n < 0, the topological recursion ensures that the W^f'^ are meromorphic 
forms with even poles only at branchpoints, and of degree at most 6g — 4 + n, we may 
decompose them on a basis of such meromorphic forms. Consider: 

5,,,„(C)= Res B{CX') (x(C')-x(a,))-"-^ 

which is a meromorphic form, whose only pole is at a^, and which behaves like: 

dz 

i?.,„(C)~(2n-l)^ + reg 
where we recall that z = a/ x{() — x{ai). 
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Therefore we can write w!f^ as a finite linear combination of Ba^ ks: 

n 

which is a finite sum since each kj is between: 

< A;^- < 3^ - 3 + n. (7-3) 

7.1 Expansion near other points 

One may also choose to expand W-!f\(i, . . . ,(n) as a Taylor or Laurent series near 
any other point p (not necessarily a branchpoint), and in powers of any other variable 
5 = /(C): 

In general, the coefficients A\^^ of Laurent expansion of a function in terms of 
one local variable Zi, are related to the coefficients A^j^^ ^ of Laurent expansion in 
terms of another variable by the Lagrange inversion formula, which just amounts 
to compute the residues: 

n 

= • • • f^es ^n\Ci. • • • , Cn) n /(^^)'"- (7-4) 

1=1 

Here, it suffices to compute the Taylor or Laurent series expansion of the basis 
forms -Bfli.n in the parameter z: 



k 



5a„n(C) = y'^a^niM ^ ^ ^ ^Z , = ReS Ba^^n{C)f{0 

k 

And that gives: 

n 

^p^!fci;...;p„,fc„ = ^ ^i!|,2mi;...;ai„,2m„ JJ^ -^ai^,mj;pj,fcj (7-5) 

ii,mi;...;i„,m„ j=l 

Remark on Cut and Join equations 

The topological recursion implies some recursive equations among the coefficients 
^a^|,fci;...;ai ^^"^ therefore, through equation Eq. (17-5^ . they imply some relation- 
ships among the coefficients A^f\. . j. . 

Those relationships can be thought of as " Cut and Join" equations. 
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Indeed, we shall see below, that for rather canonical choices of /(C)) they correspond 
to Tutte's equations for discrete surfaces, to Cut and Join equations for Hurwitz num- 
bers, or to Mirzakhani or Virasoro equations for intersection numbers of tautological 
classes. 

In fact, in most of the known applications of the topological recursion in physics and 
mathematics, some recursion relations (Cut and Join, or Tutte) are known in terms of 
the coefficients and not directly A^^^ i.e. in terms of the Laurent expansion of 

a local coordinate of the spectral curve typically near a pole or a logarithmic singularity 
(see the examples in the next section), not near branchpoints. 

7.2 Some canonical choices of expansions 
7.2.1 Meromorphic case 

It was observed in [23], that, if the form ydx is meromorphic, the coefficients of its 
Laurent series expansion near its poles, are the KP times: 

where di is the degree of the pole of ydx at pj, and ^^^ (C) is a local coordinate near pj 
given by Zp.{Q = x{() ^^'^'^^pj ^ if x has a pole at pj, and Zp^{() = x{() — x{pj) if x has 
no pole at Pj. 

This shows that a natural choice of expansion is to choose z = f{() = l/{zp.{()) 
near a pole pj. If we apply the Lagrange inversion formula as above, and expand the 
Wn^'s in Laurent series in the variables z, we shall get a natural expansion in terms of 
KP times. 

This kind of expansion is deeply related to the Frobenius manifold structure of 
moduli spaces of worldsheets. 

Example: maps, discrete surfaces, 1-matrix model 

The formal 1-matrix model provides generating functions for counting maps (maps in 
the sense of combinatorics, i.e. graphs embedded on Riemann surfaces, also called 
discrete surfaces or ribbon graphs). For the 1-matrix model, the spectral curve is 
algebraic, and x{() is a meromorphic function, with two simple poles. Moreover, the 
spectral curve is hyperelliptical, there is an involution ( ( on C, for which x{() = 
x(C), and it turns out that all stable wjf^^s are odd under that involution. This implies 
that computing the Laurent expansion at one pole of x is equivalent (up to a sign (—1)") 
to computing it at the other pole. 
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In that case we choose f{() = x{(), and we write: 

Wi^\Cu...Xn)= U"^-'--' dx,. (7-7) 

It is well known that the coefficients A^^^ are the generating functions which count 
discrete surfaces of genus with n marked faces of lengths ki, . . . ,kn and a marked 
edge on each marked face: 

#Aut(m) J- J- ^ 

where the sum is over the set of maps (or discrete surfaces) m of genus g and n marked 
faces of lengths ki, . . . ,kn and a marked edge on each marked face. nj{m) denotes the 
number of unmarked faces of m of valence j, and tj is the KP-time, found from the 
Laurent series expansion of ydx at large x: 

ydx = tj x^~^ dx. 

3 

It is also known that the topological recursion, written in terms of A^^^ , reduce 
to Tutte's equations [391 SO] for discrete surfaces (also called loop equation under their 
matrix model's representation [36] )• In particular, let us derive them for the simplest 
case of the disc amplitudes, n = 1 and (7 = 0. 

Tutte's recursion consists in removing the marked edge from the marked face of 
degree k + 1 and enumerating all possible maps obtained in this way: 

• either the other side of the marked edge is in another face of degree I and removing 
it gives a map with one marked face of degree /c + / — 1; 




• either the other side of the marked edge is the same face and removing this edge 
disconnects the map into two components: one with one marked face of degree 
/ < k and one with one marked face of degree k — I — 1. 
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This procedure is bijective and gives the following relation: 

1=0 1^2 

which can indeed be thought of as a cut-and-join relation for the discs. In terms of the 
generating function Wi^\x) = ^ + X]fc>i Af^x~^~^ dx = ydx, it reads 

y^ + yV\x)=P{x) = [y V'{x)]^ 

where V'{x) = JZj'^j^'^'^y P{x) = [/(x)]^ is the positive part of the Laurent 
expansion in x of /(x), thus P{x) is a polynomial of x. 

This implies that the disc amplitude satisfies an algebraic equation, and we have 

y=^ (V'{x) ± y/V'{xy-4P{x)') . 

This equation defines the spectral curve of this matrix model, or enumerative problem 
of maps. 

Using this procedure to remove one edge of maps of arbitrary topology gives more 
general Tutte's equations which define recursively the coefficients A^^^ ^ j it was 
proved in [20] these recursions are equivalent to the topological recursion for theiyi^^'s, 
and thus that the A^j^^ are the result of the Lagrange inversion formula on the topo- 
logical recursion for the coefficients A^^^ 

7.3 Case of non-meromorphic singularities 

If ydx has non meromorphic singularities, then f[x) =power of x can't be a good 
expansion parameter. 

For example, assume that x has a logarithmic singularity, i.e. assume that e^ has a 
meromorphic singularity at some pole p. Then, it is natural to expand in powers of e^. 

This type of logarithmic singularity occurs for applications to topological strings, 
because the spectral curve is of the form if(e'^,e^) = where if is a polynomial. In 
that case, e^ and e^ are meromorphic functions on C. 

Then, one can compute, through formula Eq. fl7-5p . the coefficients of the expansion 

in terms of the A[^^ .^^ computed at branchpoints. 

The relationship between the coefficients A^^^ which compute numbers of world- 
sheets having given perimeters near the log singularities of x, and the coefficients 
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Af: fci- -ai k„ computing worldsheets having given perimeters near branchpoints, can 
be thought of as a kind of generahzation of ELS V formula [13] . 

Also in that case, the topological recursion written for the coefficients A^^^ can 
be viewed as a generalization of the cut and join equations [301 EI]- 

Example: Hurwitz numbers and ELSV 

The spectral curve for Hurwitz numbers is related to the Lambert function y = L(e^): 

e^ = ye~y. 

This means that y is a meromorphic function on CP^, with a pole which we choose to 
be at 0: 

and X is not meromorphic, it has a pole at 2; = and a log singularity aX ( = 1 and 
C = 00: 

1 1 00 

-(0 = -l + T + ln(l-.) = -l-E^- 

^ ^ k=2 ^ 

There is a unique branchpoint a.t ( = a = 00, where dx vanishes. The local 
parameter near the branchpoint is: 

For any spectral curve with only one branchpoint, it is easy to write the expansion 
of Wi'^ near the branchpoint ( 00, i.e. 2; — * 0, in terms of intersection numbers of 
tautological classes, see [H]. Indeed, since the topological recursion computes residues 
only at the branchpoint, we need only to know the Taylor expansion of y{() — y{() 
near the branchpoint: 

yiO - yiO = 5^ W(a^(C) - x{oc))'-'^ (7-9) 

k 

and this quantity is exactly the Kontsevich spectral's curve with times h. The Wi'^'s 
can then be expressed in terms of intersection numbers of ip and n classes, see [21] . 

On the other side, it is known that the expansion near C = 1 in terms of e^ gives 
the Hurwitz numbers [TT] . 

In that case, the Lagrange inversion formula can be viewed as the ELSV formula 
[T3] , and the topological recursion in terms of Hurwitz numbers can be viewed as the 
cut and join equations of Goulden- Jackson- Vakil [30| [31]. 
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Let us point out that the spectral curve was obtained by these cut-and-join equa- 
tions for the disc amphtudes, as in the discrete surfaces case. Hurwitz numbers hg,^ 
enumerate coverings of CP^ by genus g surfaces ramified over infinity with profile fi 
and at most simply ramified anywhere else. The Riemann Hurwitz formula fixes the 
number of simple ramification points away from infinity, to be 2(7 — 2 -|- /(/i) -|- 
One can get a recursion on the Hurwitz numbers hg^^ by removing (or resolving) one 
simple ramification point from the such coverings and enumerating all possible changes 
of the ramification profile fi compatible with such a resolution. In particular, for g = 
and /i = {n}, this procedure gives cut-and-join equations defining recursively the "disc 
amplitudes" /io,n- Let us define Hg^^ := j^jjj^^^fi^hg^fj,, the cut-and-join equation 
reads: ^ 

^-f^O.n = 2 -^0,fc-f^O,n-fc- (7-10) 

^ k=l 

This can be turned into a differential equation for the generating function y = 

2y-y'= [ ^d{en (7-11) 
J X 

which is a mere rewriting of the Lambert equation y = e^+^. The cut-and-join equations 
thus allow to find the spectral curve associated to this enumerative problem in the 
patch near C = 1- The other cut-and-join equations are then obtained by the Lagrange 
inversion formula from = 1 to C = cxd and lead to the topological recursion. 

7.4 Decomposition on branchpoints 

When the spectral curve has a single branchpoint a, since the topological recursion 
computes residues only at the branchpoint, we need only to know the Taylor expansion 
of y{() — y{() near the branchpoint {z = ^y x{() — x{a)): 



y(C)-2/(C) = 2^-$^W z''-' 

k>l 



and we recognize exactly the Kontsevich spectral's curve with times = ;^ Tr A~'^ (see 
[231 [21]). The Wn^'s and Fg = W^^^'s of the Kontsevich integral can then be expressed 
in terms of intersection numbers of ip and k classes, see [21], and the result is: 



Wjf\C,, . . . Xn) = 2-'^-(2e-*°)^-2^-« Yl Et:^ E 



(7-12) 



do-\-di-\-...-\-dn=dg,n k=0 bi-\-...-\-bk=do,bi>0 

i=i 1=1 1=1 i=i 
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where Zi = \/x{(i) — x{a), and the dual times ik are related to the t^s by the following 
transformation: 

ti = , 2 - = 2e~*« 
/W = E^^^^-" - /(^) = -ln(l-/(.)) = f;t,.^ 

a=l ■ ^ 6=1 

For a given branchpoint a, we write: 

^Kontsevich(a) = 6^9-^". 

When there are several branchpoints, it was shown in [38] , following [5l [6] that the 
symplectic invariants Fg can be obtained in terms of a product of Kontsevich integrals 
at each vertex: 

= e^™ . Y[ ZKonteevich(ai), (7-13) 

ai 

where f/mixing is a mixing operator: 



with 



dzi{Ci)dzj {C2] 



and the differential operator: 

a(c):^i:«..,'(CM.,c)-^^. 

using the notations: 

Zi{0 := Vx(C) - x(a,) and t/(C) - y(C) = 2z.(C) - ^ Wi,i^f-'(C). 

fc>i 

The times tk^i are thus the Kontsevich times of the z'th factor of the expression 
Eq. (JHSD- 

This mixing formula shows that all coefficients A^^^ , , can be written in 
terms of intersection numbers of tautological classes. 

Therefore, the Lagrange inversion formula allows to express the coefficients 
A^^^ corresponding to another expansion, in terms of intersection numbers. 

This can be viewed as a generalization of the ELSV formula. 
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8 Application: topological strings and BKMP con- 
jecture 



The previous sections were aimed at general spectral curves. Here, we focus on spectral 
curves related to the mirror geometry of toric Calabi-Yau 3 folds, and the application 
to topological strings. 

In [To], Bouchard, Klemm, Marino and Pasquetti conjectured that the open and 
closed amplitudes of type A topological string theories on some Toric Calabi-Yau 3- 
folds coincide with the symplectic invariants of their mirror B-model's target space. In 
this section, we explain why the computation of the Gromov-Witten invariants (closed 
and open - or relative as defined by [35]) reduces to the enumerative problem solved 
earlier in this paper. Once again, this is not a proof but just a heuristic explanation 
of this conjecture. 

8.1 Toric Calabi-Yau 3- folds and localization 

Let us consider a A-model topological string theory whose target space X is a toric 
Calabi-Yau three-fold. 

One of the main features of this type of geometry is that it can be realized as the 
gluing of a set of patches, as a fibration over a non-compact convex subspace 
of M^. Its geometry can be described by a Toric graph together with the Kahler 
parameters of X. The Toric graph corresponds to the degeneration of the fibers 
(of the fibration): the lines of the Toric graph represent the locus where two out 
of the three 5'^'s shrink to zero (see figHJ] for the simplest example X = C^). In 
particular, the external legs of the Toric graph correspond to such degeneration along 
a non-compact direction: such external legs are thus C's which can be compactified to 
CP^'s by including the point at infinity. 

This target space can be equipped with many possible Torus actions. Each of them 
leads to a different localization computation and thus, a different hamiltonian system. 
In order to match the classical computations of topological string theories, let us choose 
one particular S'^-action on X. For this purpose, choose one external leg L23 of the 
Toric graph and denote by X2 and X3 the two coordinate of X which vanish along this 
edge (they correspond to the two shrinking S'^'s along this external leg). Denote by 
Xi the remaining coordinate of X which vanishes along the two other legs L12 and L13 
of the vertex at which L23 ends. 

Now consider the Torus action on C^: 

(Xi,X2,X3) ^ (t''^Xi,t^^X2,t''^X3). 
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One of the fixed points of X under this Torus action is the tip of L23 defined by 
Pi := {Xi ^ 00, X2 = 0,X3 = 0). If pi = f, P2 = —f and ps = for some integer /, 
this action reduces to a S^-action which we shall consider from now on. 




Figure 15: To the left: Toric graph of C^. The three planes Xj = intersect along the 
edges of the toric graph. The black circle represents a level of the Torus action with 
Pi = /, P2 = — / and ps = for some integer /: {\Xi\ = r, X2 = 0,^3 = 0) for some 
constant r. To the right: the same graph with the tip pi of the external leg included 
to compactify the corresponding direction. 



Let us now enumerate holomorphic maps: 



'^g,n ■ ^ X 

from a genus g Riemann surface with n marked points to X, i.e. compute integrals 
of the form 

Aid) - 



where M.g,n is the compactification of the moduli space of such maps (see [33] for 
example). 

The S'^-action on X induces a hamiltonian action on the moduli space Mg,n- The 
Atiyah-Bott localization formula states that such integrals reduce to integrals over the 
fixed locus of the Torus action: 



where {Fixg^n} runs over the fixed locus of the Torus action and e{Np^^^ n/Mg J 
Euler class of the normal bundle of Fix„ „ in 



rigorous approach requires the introduction of an obstruction theory and consider the integrals 
over the virtual fundamental class built out of it. Nevertheless, the fixed locus is the same as for the 
compactified moduli space: only the integrant is changed in the localization formula (see [32]). 
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This means that the observables reduce to a sum over all maps whose images in X 
are invariant under the hamiltonian S^-action counted with a weight ^{Np^^^ n/Mg 
In particular, the marked points have to be mapped to fixed points of the S^-action in 
X by maps in FiXg^n- 

For a fixed point € FiXg^n, '^~g\i allows to lift the orbits of the Torus action in 
X to closed paths generating S^: the level lines of the S'^-action equip the worldsheet 
Tjg with a complex structure. The coordinates in this foliation of the worldsheet are 
nothing but the action-angle variables of the Hamiltonian system under study, the 
circles are the horizontal trajectories of this hamiltonian action, at least locally near 
the point pi := {Xi = oo, X2 = 0, X3 = 0). 

One has thus reduced the computation of A^f^'s to the study of embeddings of world- 
sheets into a Jacobian thanks to action-angle variables. This is exactly the problem 
studied in the previous section. 

One can refine the description of the fixed locus of the action in M.g^n- For such 
a Torus action, one can decompose An^ by fixing the image of the marked points as 
well as the local behavior of the enumerated maps. Let pi, . . . ,Pm be fixed points of 
the S'-'^-action in X. For a set of partitions /ii, . . . ,/im, let -FzXg^^^^^j,...,^^ be the set of 
maps mg,^i,^2,.--,Mm such that Z(/Xj) marked points are mapped to pi 

1=1 

with ramification profile /Xj. We finally denote 

the localized integrals restricted to this subset of the fixed locus. 

Remark that these amplitudes are closely related to "open" amplitudes which are 
integrals over the moduli space of maps from open worldsheets to X with fixed Brane 
boundary conditions by removing a small circle from the closed surface around the 
marked points. 

Let Eg be a genus g surface with n marked points zi, . . . ,Zn embedded into X thank 
to a stable map mgn G FiXg^^ with /(/i) = n and Zi, . . . ,Zn all mapped to the fixed 
point pi := {Xi = 00, X2 = 0, X3 = 0) of the compactification of X. 

For r G M large enough, the pull-back of the circle {\Xi\ = r, X2 = 0,X3 = 0) 
by m~l^ has n connected components winding /ij times around Zi respectively. Thus, 
removing the pull-back of the discs {\Xi\ > r, X2 = 0,^3 = 0) by from 
leaves us with an open surface E^ „ with n boundaries whose embedding in X by rrig n 
is stable under the Torus action with its boundaries mapped to the Brane {\Xi\ = 

^"This localization formula is the sum over instantons of sec l2.3l One instanton is a map stable 
under the Torus action, i.e. a fixed point of the moduli space. 
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Figure 16: Example of stable map contributing to ^23. The blue (resp. red) cycle on 
the worldsheet winds twice (resp. three time) around the fixed point pi in the target 
space. 



r,X2 = 0,^3 = 0). This mapping between stable open surfaces and marked closed 
ones is bijective up to symmetry factor and allows to get 



A(f)ocAr(^)(r):= Yl I 



where Fixg^^{r) stands for the fixed locus of the moduli space of maps from an open 
surface with /(/i) boundaries mapped to the brane {\Xi\ = r, X2 = 0, X3 = 0). N'l^\r) 
are exactly the volume computed in the preceding sections. 

Remark 8.1 One can invert this procedure by recapping the open surface by a disc (or 
rather a semi- infinite cylinder). This gives the precise relation between the open and closed 
amplitudes. This is the relation I5-16[ 

It is also interesting to remark that in this picture, the length of the boundaries of 
the open surfaces are quantized once r is fixed: they are labeled by the winding num- 
bers. The Laplace transform considered for the definition of the generating functions 
of these numbers in section [7] are thus replaced by discrete Laplace transforms with 
respect to the winding number J"]. 

8.2 Cut-and-join, spectral curve and mirror map 

Following Kontsevich [33] , one can describe the fixed locus of the moduli space of maps 
FiXg^fj,-^^^^^,,,^^^ rather explicitly by associating a decorated graph to each component of 
this fixed locus. In particular, the case 5^ = 0, /^i = {ci} and /ij = for i ^ 1 maps to the 



These transforms are the one used in [Tni[IH[lHll3] for the computation of Hurwitz numbers and 
Gromov-Witten invariants of C^. 
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enumeration of rooted trees with labeled edges. The integrals of the Euler classes on 
each component can then be performed explicitly (see [15] for the framed vertex or [32] 
for the general Toric CY 3-fold case) and one ends up with a problem of enumeration 
of trees with fixed weight. 

This procedure gives the disc amplitudes. The generating function for these disc 
amplitudes was proved (at least physically [1]) to be closely related to the superpotential 
of the mirror theory: the problem of enumeration of trees can be solved by induction 
on the winding number the generating function Wi'^\x) for these numbers is thus 
a "tree function" solution to an equation of the type if(e^,e^i ''^^^) = which is the 
equation of the singular locus of the B-model target space. 

Since the disc amplitudes are the only unknowns of the topological recursion, once 
we know the disc amplitudes, all the other amplitudes can be computed by the topo- 
logical recursion. As explained in the preceding section, the generating function for the 
disc amplitudes defines the spectral curve. One thus recovers that the Gromov-Witten 
invariants of some A-model are given by the symplectic invariants of the corresponding 
B-model target space in the coordinates associated to A-Branes under study. 

Many questions arise from this approach and have to be addressed in order to clarify 
the combinatorics underlying the topological recursion. An important one is to make 
clear how the choice of a torus action maps to the choice of parameterization of the 
B-model target space (or of the spectral curve in our language). 

Given a Torus action used for the localization method, there is no ambiguity in 
the determination of the integrals over the moduli space of maps: it fixes both the 
polarization and the framing by fixing the Torus action. In the simplest case of a C'^ 
target space, the map between the choice of a Torus action and a parametrization of 
the spectral curve has already been performed (see for example [15] or [12]) and the 
framings of the corresponding Toric vertex. 

For a general toric CY target space, the theory of the topological vertex [53] allows 
us to follow the same procedure and associate a parameterized curve to a given torus 
action, i.e. it maps a foliation of the worldsheet to a parameterization of the spectral 
curve. The topological recursion follows from a particular local parameterization of 
the spectral curve. It would be interesting to see if this local parametrization of the 
spectral curve can be mapped to a local action on the target space, and thus, to a 
local foliation of the worldsheets. 
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9 Conclusion 



First, let us say again that all the construction presented in this paper is semi- heuristic 
and many details need to be further studied and made precise in order to become 
really rigorous. It can't be seen as a proof of the "remodeling the B-model" formalism, 
but rather as an intuitive geometric understanding of it. One could expect to find a 
rigorous proof of BKMP by using generalized cut and join equations, as we suggest 
in section [TJ this is more or less what was done for Hurwitz numbers in [19] and for 
X = by Zhou [H US] and Chen [H]. 

Reverse engineering 

Our approach is a reverse engineering of the topological recursion. In fact, knowing 
that a string theory satisfies the topological recursion, and knowing that the topological 
recursion depends only on the data of a spectral curve, we tried to reconstruct the string 
theory from the spectral curve. The topological recursion clearly implies that every 
worldsheet should be decomposable in a unique way into propagators and cylinders. 
Therefore, there is no really other choice than the theory we presented here. 

The main drawback of that reverse engineering approach, is that, given a string 
theory, it is not so straightforward to recover the spectral curve associated to it. 

In the context of topological B-strings in toric CY 3-fold target spaces, the spectral 
curve was found from mirror symmetry [32], and the "remodeling the B-model" idea 
of BKMP [TO] is based on that. 

For other models, it would be interesting to see how our method applies in practice. 
For instance, starting from the Lambert spectral curve, how do we see that our foliated 
flat worldsheets are realizing branched coverings of the sphere. 

Also for matrix models, the spectral curve is well known, it is algebraic, it is re- 
lated to the Toda chain integrable system, and it is known that it enumerates discrete 
surfaces (called "maps" by combinatorists). It would be interesting to see how our 
fiat coordinates foliate discrete surfaces in that context. In other words, what are the 
horizontal and vertical trajectories on the discrete surfaces ? 

As pointed out in the preceding section, for the works on Hurwitz coverings and 
matrix models, the spectral curve seems to always come from a cut-and-join procedure 
(or Tutte's equations for the matrix models) recursively defining the disc amplitudes. 
The latter takes the form of an induction of a set of trees whose weights depend on 
the propagator amplitude, i.e. the possible singularities of the chosen hamiltonian 
fibration. The study case by case of this general statement would be very interesting. 
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Flat coordinates and flat surfaces 

In fact, this idea of foliating worldsheets with the help of the fiat connection of an 
integrable system was already used in string theory. For instance in [T7|, Dorey et 
al. used a fiat connection to parametrize bare cylinders. They didn't consider higher 
topologies. 

Let us also mention the linli with the theory of fiat surfaces developed by Zorich et 
al. It is clear that our fiat surfaces, can be realized as fundamental polygonal domains 
of a plane modulo the lattice, l.G. clS db polygonal (embedded in the Jacobian, see 
figJTTI). with some opposite sides glued together. There is an important literature on 
that theory of fiat surfaces, see [H], and it would be interesting to clarify the linl<; with 
our approach. 





Figure 17: The linear constant motion in the periodic Jacobian generates a fiat surface 
obtained by gluing opposite sides of a polygon. The boundary at time t = has 
the topology of a circle, but at some critical time, it ceases to be a circle, there is 
a branching. The image of this constant motion, produces a Riemann surface in the 
target space X. 



Integrable systems 

The theory of Hitchin's systems, relies on a foliation of moduli spaces, with the orbits 
of an integrable system. The foliation concerns the moduli space itself, and not the 
worldsheets, however, it is clear that there is a link between those two approaches ( a 
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foliation of worldsheets induces a foliation of moduli spaces), and it would be worth 
developing it. From this perspective, it would be interesting to compare our approach 
with the description of topological string theories by Gerasimov and Shatashvili in [29] . 

Another link between enumerative geometry problems and integrable systems arises 
through the Frobenius manifolds structure. There also, it seems interesting to under- 
stand how our approach fits in the framework of Frobenius manifolds. In particular, 
the relation between flat and canonical coordinates in the work of Dubrovin and al. 
[TS] seems to coincide with the transform expressing the expansion in powers of local 
coordinates near the branch points, in terms of the expansions in terms of KP times 
related to the expansion near the poles of ydx, see section [71 This points towards the 
link between a change of torus action in the topological string setup and the structure 
of Frobenius manifolds. 

Other prospects 

In [16], Dijkgraaf and Vafa showed that Koddaira-Spencer theory also satisfies the 
topological recursion, and thus there is a quantum field theory equivalent to this for- 
malism. Another way of seeing a quantum field theory, is through integrability. In- 
tegrable systems' correlation functions are determinants, and can be written as free 
fermions integrals. In particular, it would be interesting to generalize Kostov's work 
on CFT description of matrix models [M] to realize the correlation functions studied 
in this paper as correlation functions of some associated conformal field theory on a 
Riemann surface with insertion of Twist operators associated to the branch points. 

Finally, we considered only classical integrable systems in these notes. In [T3], the 
topological recursion has been generalized to quantum spectral curves (C, x, y) with 
non-commuting x and y: [y, x] = h. It seems very likely that this quantization is related 
to Nekrasov's partition function for equivariant theories [37] . Generalizing the method 
developed here in this non-commutative context could give a new geometrical point 
of view on this partition function as well as cut-and-join like formulas for equivariant 
theories with higher dimensional torus actions. This could also lead to a link with the 
AGT conjecture [1]. 

Also, let us mention that here we considered only "closed" strings, ending on D- 
branes, with Dirichlet boundary conditions, i.e. horizontal trajectories. It seems easy 
to extend our method to worldsheets bounded also by vertical trajectories, which thus 
correspond to Von Neumann boundary conditions, and which realize open strings. It 
seems that one could easily extract a topological recursion formula for open strings as 
well. 
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